Resonating valence bond states in the PEPS formalism 
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We study resonating valence bond (RVB) states in the Projected Entangled Pair States (PEPS) 
formalism. Based on symmetries in the PEPS description, we establish relations between the toric 
code state, the orthogonal dimer state, and the SU(2) singlet RVB state on the kagome lattice: We 
prove the equivalence of toric code and dimer state, and devise an interpolation between the dimer 
state and the RVB state. This interpolation corresponds to a continuous path in Hamiltonian space, 
proving that the RVB state is the four-fold degenerate ground state of a local Hamiltonian on the 
(finite) kagome lattice. We investigate this interpolation using numerical PEPS methods, studying 
the decay of correlation functions, the change of overlap, and the entanglement spectrum, none of 
which exhibits signs of a phase transition. 



I. INTRODUCTION 

Resonating valence bond (RVB) states have been intro- 
duced by Anderson^ as a wavefunction in the context of 
high-temperature superconductivity and have since then 
received significant attention. RVB wavefunctions are 
obtained as the superposition of all nearest-neighbor (or 
otherwise constrained) singlet coverings on a given lat- 
tice, and they have been studied as an ansatz captur- 
ing the behavior of frustrated quantum spin systems, in 
particular Heisenberg antiferromagnets on frustrated lat- 
tices.— They are believed to describe so-called topological 
spin liquids, i.e., exotic phases with degenerate ground 
states which however do not break rotational or transla- 
tional symmetry, despite the presence of strong antiferro- 
magnetic interactions. Particular interest, both theoreti- 
cally and experimentally, has been devoted to frustrated 
magnets on the kagome lattice, as those systems are real- 
ized in actual materials and form candidates for the first 
experimental observation of a topological spin liquid.— — 

Unfortunately, models with potential RVB ground 
states are typically difficult to study. One the one hand, 
this is due to presence of frustration. Yet, the bigger issue 
seems to be that different singlet configurations are not 
orthogonal, which up to now has e.g. hindered a full un- 
derstanding of how RVB states appear as ground states 
of local Hamiltonians.^ In order to better understand 
the RVB phase, simplified so-called dimer models have 
been introduced^ where the system is re-defined on the 
Hilbert space spanned by all dimer coverings of the lat- 
tice, making different dimer configurations orthogonal by 
definition. Dimer models have subsequently been studied 
for various lattices, and it has in particular been found 
that for certain geometries, such as for the triangular ^^ or 
the kagome lattice, ^^ such states appear as ground states 
of local Hamiltonians with topological ground state de- 
generacy (four- fold on the torus), forming a Z2 topolog- 
ical spin liquid. Unfortunately, these findings cannot be 
mapped back to the true RVB state due to the differ- 
ent underlying Hilbert spaces which are not related by a 



simple mapping, and thus, the way in which the corre- 
sponding RVB states can describe ground states of local 
Hamiltonians is only partially understood (although very 
remarkable progress in this direction has been made dur- 
ing the last year^ii^). 

More recently. Projected Entangled Pair States 
(PEPS) have been introduced as a tool to study quan- 
tum many-body wavefunctions . ^^1^^ PEPS give a descrip- 
tion of quantum many-body states based on their entan- 
glement structure in terms of local tensors, and form a 
framework which allows to both analytically understand 
these wavefuctions, and to study their properties numer- 
ically. In particular, there is a clear way to understand 
how PEPS arise as ground states of local Hamiltonians, 
especially for systems with unique ground states^ ^ and 
with topological order,— based on the symmetry prop- 
erties of the underlying tensor representation. Numer- 
ically, they form a tool for computing quantities such 
as correlation functions or overlaps efficiently with very 
high accuracy, using transfer operators of matrix product 
form.— 

In this paper, we perform a systematic study of RVB 
and dimer states in the PEPS formalism, focusing on the 
kagome lattice. We introduce closely related PEPS repre- 
sentations for the RVB wavefunction^^ and a version of 
the dimer state in which different dimer configurations 
are locally orthogonal. These representations allow us 
to derive a reversible local mapping between the RVB 
and the dimer state, as well as to prove local unitary 
equivalence of the dimer state and Kitaev's toric code.— 
This yields a Hamiltonian for our dimer model, and sub- 
sequently a Hamiltonian for the RVB state: This is, we 
can for the first time prove that the RVB state on the 
kagome lattice is the ground state of a local Hamiltonian 
with topological ground state degeneracy (four-fold on 
the torus) for any finite lattice. The PEPS formalism 
further allows us to construct an interpolation between 
RVB and dimer state where the Hamiltonian changes 
smoothly along the path. Using techniques for numer- 
ical PEPS calculations^ we try to assess whether the 




FIG. 1. a) Oriented kagome lattice. Spins are associated 
to vertices, b) Dimer covering of the kagome lattice. The 
lattice is completely covered with dimers (marked blue), i.e., 
disjoint pairs of adjacent vertices. 



RVB state is in the same phase as the dimer state (which 
is equivalent to the toric code and thus a Z2 topological 
spin liquid) by looking for signatures of a phase transi- 
tion along the dimer-RVB interpolation. We have con- 
sidered the behavior of correlation functions, the rate at 
which the ground state changes in terms of the wave func- 
tion overlap,^^ and the entanglement spectrum of the sys- 
tem, and have found that all of these quantities behave 
smoothly and show no sign of a phase transition. 

Let us describe the organization of the paper. For the 
sake of conciseness, alternative definitions, proofs, etc. 
have been moved to appendices. In Sec. [Ill we introduce 
the RVB and dimer states and their PEPS representa- 
tions, as well as the formulation using tensor networks. 
Alternative PEPS representations are discussed in Ap- 
pendix El In Sec. mil we discuss the relations between 
the toric code, the dimer state, and the RVB state, us- 
ing the symmetry of the underlying tensors. Based on 
these findings, we show how to construct a smooth in- 
terpolation between the dimer state and the RVB state 
in Sec. IIVI and subsequently study its properties numer- 
ically. In Sec. |Vl we use the relation between the toric 
code, the dimer state, and the RVB state to construct 
parent Hamiltonians for the dimer state and the RVB, 
as well as a smooth path of Hamiltonians for the inter- 
polation between them. Appendix O provides a simpler 
parent Hamiltonian for the dimer state, based on a more 
direct mapping to the toric code, and Appendix [Dl shows 
how to directly derive a (much more compact) parent 
Hamiltonian for the RVB state by generalizing the tech- 
niques developed in Ref. [l6| for PEPS with symmetries 
to the symmetries found in the RVB. 



II. DEFINITIONS 

A. The RVB and orthogonal RVB state 

We start by introducing dimer states and resonating 
valence bond (RVB) states. We focus on the kagome 




FIG. 2. PEPS construction of the RVB state: Place |e) states 
(green), Eq. ([T]), and then apply the map V (red), Eq. (|2|, 
as indicated. To obtain the dimer state, replace V by Pi, 
Eq. da)). 



lattice. Fig. [TJi, both for its relevance and for clarity of 
the presentation; but our techniques generalize to other 
lattices (see conclusions). A dimer is a pair of vertices 
connected by an edge. A dimer covering is a complete 
covering of the lattice with dimers. Fig. [Ud. We can as- 
sociate orthogonal quantum states |D) with each dimer 
covering D. Then, the dimer state is given by the equal 
weight superposition l^dimer) = Y.\D)^ where the sum 
runs over all dimer coverings D. Usually, it is favorable 
to ensure orthogonality of different states \D) and \D') lo- 
cally; we will introduce such a version of the dimer state 
which is particularly suited for our purposes in Sec. Ill CI 
Let us now turn towards the resonating valence bond 
(RVB) state. We first associate to each vertex of the 
lattice a spin-| particle, or qubit, with basis states |0) = 
It) and |1) = |i). Then, for each dimer covering D we 
define a state \(j{D)) which is a tensor product of singlets 
|01) - 1 10) (we omit normalization throughout) between 
the pairs of spins in each dimer in the covering, where the 
singlets are oriented according to the arrows in Fig. [TJi. 
The resonating valence bond (RVB) state is then defined 
as the equal weight superposition |^rvb} = Y.dW{^)) 
over all dimer coverings. 



B. PEPS representation of the RVB state 

We will now give a description of the RVB state 
in terms of Projected Entangled Pair States (PEPS).— 
PEPS are states which can be described by first placing 
"virtual" entangled states between the sites of the sys- 
tem, and subsequently applying linear maps at each site 
to obtain the physical system. A PEPS representation of 
RVB states has first been given in Ref. [13; a detailed dis- 
cussion how it is related to our description can be found 
in Appendix [Al 

To obtain a PEPS description of the RVB state, we 
first place 3-qutrit states 



k)= E eijk\ijk) + \222) 

i,j,k=0 



(1) 



inside each triangle of the kagome lattice, as depicted in 



Fig. [21 Here, Sijk is the completely antisymmetric tensor 
with £012 = I7 and i, j, and k are oriented clockwise (i.e., 
consistent with the arrows in Fig. [1^). Second, we apply 
the map 



P = |0)((02| + (20|) + |l)((12| + (21|) 



(2) 



at each vertex, which maps the two qutrits from the ad- 
jacent l^} states to one qubit. It is straightforward to 
check that this construction exactly gives the resonating 
valence bond state defined above (see Appendix [Al for a 
detailed discussion). 



C. PEPS representation of dimer state 

In a similar way as for the RVB state, we can also 
obtain a PEPS representation of the dimer state. To this 
end, we enlarge our local Hilbert space and replace the 
map V in Fig. [2] by the map 

V^ = |02)(02| + |12)(12| + |20)(20| + |21)(21| . (3) 

It is straightforward to check that the resulting state is 
an equal weight superposition of all dimer coverings, and 
that different dimer configurations are orthogonal, since 
the position of the dimers can be unambiguously inferred 
from the location of the |2}'s. Note that while we write 
the physical space as a subspace of a 3 x 3-dimensional 
space, in fact only a four-dimensional subspace is used. 
A more detailed discussion of the construction can again 
be found in Appendix [Al 

As one might expect, it is possible to obtain the RVB 
state from the dimer state by coherently discarding the 
information about where the singlet is located and only 
keeping the singlet subspace span{|0}, |1}}. This is most 
easily seen by noting that 

P = VV^ , 

and thus, applying V to each site of the dimer state yields 
the RVB state. 

Note that our representation of the dimer state differs 
from the formulations typically used in the literature in 
two points: Firstly, while dimer models have been usually 
studied on the space spanned by valid dimer configura- 
tions, our formulation allows to locally check whether a 
dimer is present on a given link. This allows us to en- 
force the restriction to the space of valid dimer configu- 
rations using a local Hamiltonian, and ultimately enables 
us to map our dimer model to the toric code by local uni- 
taries. Secondly, unlike mappings where the presence or 
absence of a dimer is represented by a spin,^^ our map- 
ping uses actual singlets built of two entangled particles. 
This allows us to locally interpolate between the dimer 
and the RVB model, and gives rise to a gapped entan- 
glement spectrum, which is in contrast to the gapless en- 
tanglement spectra which had been observed previously 
for topologically ordered systems. 



D. Gauge transformation 

The PEPS representation of the RVB and dimer state 
are built using unnormalized singlets, as the weight of 
|01} - 1 10} and |22} in ([1]) is the same. Since the number 
of singlets in all dimer configurations is the same, this is 
not an issue; however, in some situations (cf. Sec. IIIIBp 
it will be more convenient to normalize the singlets, i.e., 
replace ([1]) by 



^) = ^ E e,,k\ijk) + \222) 



(4) 



i,j,k=0 



Note that \i) = (F F Y)\e), where Y = 
diag(2-^/'^,2-^/'^,l). Since on the other hand, P^(Y <S) 
Y) = 2~^'^V]_ (and thus also for V^ as well as the interpo- 
lation defined in Sec. lIVp . we find that these two PEPS 
do not only represent the same total state, but are in 
fact related by a local gauge transformation, allowing to 
directly relate all properties of their PEPS description. 



E. The toric code 

Consider a square lattice of qubits, with an A-B- 

pattern on the plaquettes. Fig. [3^. Kitaev's toric code 

as 



state is the ground state of the Hamiltonian 
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(5) 



where the sums run over all A and B type plaquettes pa 
and pb^ respectively, and the tensor products of Paulis 
act on the qubits adjacent to each plaquette. 

The toric code state can be written as a PEPS^^ with 
bond dimension D = 2. A particularly useful represen- 
tation is obtained by blocking the qubits in 2 x 2 blocks 
around A type plaquettes, as indicated in Fig. [3^ and 
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FIG. 3. Construction of the toric code as a tensor network, 
a) The toric code is defined on a square lattice with two types 
of plaquettes. By taking 2x2 blocks over A-type plaquettes, 
a particularly convenient tensor network description can be 
found, b) Tensor for the toric code; depicted are the non- 
zero entries. The four auxiliary indices are in the X basis, 
while the physical sites are in the Z basis. 
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FIG. 4. Tensor network notation for the RVB state, a) Ten- 
sors are denoted as boxes with one leg per index. Connecting 
legs corresponds to summing over that index, panel a) e.g. 
shows XI /3 ^L/3^^^- b) Tensor network description of the RVB 
state; the triangluar link is shorthand for the tensor E implic- 
itly defined in Eq. (|8]), see panel c). Depending on whether 
we choose P as in Eq. (|9]) or replace it by Pi, Eq. (p^ . we 
obtain a tensor network for either the RVB or the dimer state; 
in the latter case, the physical indices ifc form double indices 



letting the PEPS projector b< 



.16 



a, /3, 7, 5=0 

for each of these blocks. Here, the sums are modulo 2, 
and the ordering of the indices is illustrated in Fig. [Sb; 
the virtual qubits are expressed in the X eigenbasis |0} = 
(|0) + |l»/x/2, |i) = (|0)-|l))x/2. 



F. Formulation using tensor networks 

PEPS can also be expressed in terms of tensor net- 
works; while both descriptions are mathematically equiv- 
alent we will utilize either of them when more convenient. 

A tensor network state is a state 






In, 



5^iv) 



(7) 



where ci 



can be expressed efficiently by a tensor net- 



work, i.e., as the sum over ai, . . . , aM of product of ten- 
sors with indices n, . . . , ^at and ai, . . . , aM^ where each 
of the tensors only has a few indices. Tensor networks 
are often expressed graphically, with tensors denoted by 
boxes with legs corresponding to the indices, where con- 
necting legs corresponds to contracting over the corre- 
sponding index, cf. Fig. |4^. 

It is straightforward to see that the simplified PEPS 
representation of the RVB state can be expressed using 
two types tensors, E and P, defined such that 
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FIG. 5. Z2-invariant tensor network, a) A tensor is called Z2- 
invariant if it is invariant unter some Z applied to all virtual 
indices, where {1,Z} forms a representation of Z2. b) Z2- 
invariance is stable under concatenation of tensors. 



which form the tensor network of Fig. |4]3 (where we use 
the shorthand Fig. St for E); it has open indices ii, . . . , ^at 
and yields the coefficient Q^^...^^^ in Eq. ([7j). Again, we 
can replace P by P^ defined by 

^x= Z(PrY;lJid)iv,^\ (10) 

to obtain a PEPS representation of the dimer state. 



III. RELATIONS 

In this section, we will show how one can transform 
reversibly between the toric code, the dimer state, and 
the RVB state. To this end, we will start by recalling 
the concept of G-injectivity of PEPS^ and subsequently 
apply it it to prove equivalence of the toric code and the 
dimer state, and to devise a mapping which transforms 
between the dimer and the RVB state. 



A. Z2— injectivity 

We start by introducing Z2 -injectivity, which has been 
introduced (as G-injectivity, with G any finite group) in 
Ref. [l6| as a tool to characterize PEPS with topological 
properties. Let Z be a unitary such that {t,Z} forms a 
representation of Z2, i.e., Z'^ = t. A Z2-invariant ten- 
sor network consists of tensors which are invariant under 
the symmetry action on all virtual indices simulateously, 
as depicted in Fig. [5^ (we restrict to equal representa- 
tions on all indices). This property is stable under block- 
ing, i.e., two (or more) Z2-invariant tensors together are 
still Z2 -invariant, see Fig. [5]3. A Z2 -invariant tensor is 
called Z2-injective if the corresponding PEPS projector 
V is injective on the invariant subspace (the subspace 
belonging to the trivial representation), i.e., if it has a 
left inverse V~^ such that V~^V is the projection onto 
the invariant subspace. Differently speaking, this means 
that there exists a linear map which we can apply to 
the physical system to obtain direct access to the virtual 
system, up to a projection onto the symmetric subspace 
n^. Z2-injectivity is stable under concatenation: When 
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FIG. 6. Transformation of the dimer state into the toric code, 
a) Pi can be rewritten as a tensor network with bond dimen- 
sion 2. b) Network obtained by replacing Pi in Fig. \^ by 
the tensor network of panel a). The tensor network can be 
grouped into triangular tensors which contain only A or B 
tensors, and which are connected by bonds of dimension 2. 



composing two Z2-injective tensors, we obtain another 
Z2-injective tensor^ 

The key point to note is that any two Z2 -inject ive 
tensor networks with identical geometry and symmetry 
representation, but different tensors, are locally equiva- 
lent: There exists a map acting on the physical system 
of individual tensors which transforms the networks into 
another. This follows since both PEPS projectors, Va 
and Vb^ have left inverses which yield the projector E^ 
onto the symmetric subspace, V^Va = V^Vb = H^, and 
thus, Va = {VaV^)Vb-, and vice versa. 

An important special case are 1^2 -'isometric tensors 
where 'P is a partial isometry, V~^ = V^: Any two Z2- 
isometric tensors Va and Vb can be transformed into 
each other by a unitary acting on the physical system. 
Again, Z2-isometry is a property which is stable under 
concatenation.— 



B. Toric code and dimer state 

In the following, we will show that the Toric Code state 
and the orthogonal dimer state are both built of Z2- 
isometric tensors, which implies that they can be trans- 
formed into each other by local unitaries. (This relation 
has been observed previously by mapping the presence or 
absence of dimers to a spin degree of freedom. ^^) In this 
section, we will consider the Z2 representation {1,Z}, 
with Z = (i_i). 

For the toric code tensor Eq. (J6j), it is straightforward 
to check that it is invariant under Z®^, and acts isomet- 
rically on the invariant subspace, i.e., is Z2-isometric. In 
order to see the same for the dimer state, we first rewrite 




FIG. 7. Non-zero orthogonal virtual states and the corre- 
sponding orthogonal physical states for the dimer state ten- 
sors (the red lines denotes singlets). The A (B) type tensors 
are supported on the odd (even) parity subspace. 



Pi, Eq. (H 



as 



Pi 



k=0 



(11) 



where Aq = Bi = |0)(0| + |1)(1|, and Ai = Bo = |2)(2|. 
Using this decomposition, we replace each tensor P^ in 
Fig. IDd with a pair of tensors A and B connected by a 
two-dimensional bond, cf. Fig. [6^, oriented as in Fig. [6)3. 
This orientation allows to block the tensors into trian- 
gles, labelled A and 5, containing only A and B type 
tensors, cf. Fig. [GJd; the triangles are connected by bonds 
of dimension two. 

The bonds are associated to vertices; a bond state |0} 
(|1}) indicates that the dimer lies inside the A (B) trian- 
gle (similar to the "arrow representation" of Ref . \2u) . As 
each triangle contains either no or one dimer, the A (B) 
tensors are odd (even) under Z®^ symmetry, and more- 
over isometric on the (anti) symmetric subspace if we use 
the gauge of Eq. (|3]); the possible mappings from bond 
to physical configurations are given in Fig. [71 

We now group one A and one B tensor to obtain square 
tensors, which are antisymmetric and isometric; grouping 
two of these square tensors finally yields a Z2-isometric 
tensor with symmetry representations Z and Z (8) Z, re- 
spectively. Blocking two tensors of the toric code also 
yields a Z2-isometric tensor with the very same symme- 
try, and thus, the dimer state and the toric code can be 
converted into one another by local unitaries. 



C. RVB and dimer state 

Let us now show that we can reversibly (though not 
unitarily) transform the dimer state to the RVB state. 
In principle, the forward direction is obvious, as 7^ = 
VV^] however, this transformation cannot be inverted. 
In order to construct an invertible transformation, we 
will show that both states can be understood as being 
constructed from Z2-injective tensors connected by 1^) 




b) 



• \^ ^M^ • • \^ ^Jt__ 



. X J 




FIG. 8. Definition of the tensor P^. a) P^ is obtained by 
blocking 12 tensors P with 6 triangular "bond tensors" E; 
it maps 6 virtual qutrits to 12 physical qubits. b) Effective 
depiction of the tensor P^ . 



bonds: Following the reasoning of Sec. IIII A\ we can re- 
versibly transform between the two Z2-injective tensors 
and thus between the RVB state and the dimer state. 

For this section, we will use the PEPS representation 
Fig. IHof both the RVB and the dimer state. The repre- 
sentation of the Z2 symmetry is given by 



4 



Z:= 



-ll 



It is straightforward to see that both P, P^, and the 
bond tensor E are anti-symmetric under the action of 
the symmetry on all virtual indices. Let us now form 
a new tensor by blocking all the tensors covering one 
star, cf. Fig. [51 We call the new tensors P^ and P^ ^, 
respectively (with corresponding PEPS projector V^ and 
T^i,^)' Clearly, we can tile the lattice with P^ tensors 
connected by E bonds, see Fig. [9l Due to stability of Z2- 
invariance, one finds that P^ and P^ ^ are Z2 -invariant. 
On the other hand, it can be checked analytically using 
a computer algebra system that V^ is indeed invertible 
on the invariant subspace of Z®^, i.e., V^ is Z2-injective. 
This immediately implies Z2-injectivity of P^ ^ (as V^ = 

We can now use the argument of Sec. IIII Al to con- 
struct an invertible mapping: Z2 -inject ivity of V^ im- 
plies the existence of V^ such that V^V^ is the pro- 
jector onto the Z® ^-invariant subspace on the virtual 
qutrits along the boundary of the star. If we now de- 
fine O^ = V^^^V^, it follows that O^V^ = V^^^, i.e., it 
(locally) converts the RVB into the dimer state. [Note 
that O^ is only well-defined on the range of 'P^, rg('P^); 
we choose O^ to vanish on its orthogonal complement.] 
Analogously, we can define O^ such that O^V^^^ - V^ 
[again such that it vanishes on rg('P^ o)^]' ^^^ ^^ ^^ easily 
verified that O^O^ = IrgCP^) and O^O^ = IrgCP^^), 
i.e., O^ is the pseudoinverse of O^. 




FIG. 9. Complete covering of the kagome lattice with stars, 
Fig. [HI The stars tensor P^ (or projectors V^), together 
with triangular tensors E (or l^) bonds), describe the RVB 
state (or its orthogonal version, if we choose P^ ^ instead). 
By acting with O^ on each of the stars, we can reversibly 
convert between the RVB and the orthogonal RVB. 



Applying the maps O^ or its inverse to all K = N/12 
stars in the tiling Fig. [9l we can now reversibly con- 
vert between the dimer and the RVB state, l^dimer) = 
C^i^l^RVB) and I^rvb) = (O^i)^^l^n 



^oRVB/ 



IV. INTERPOLATION 

A. Interpolation between dimer model and RVB 

Let us now show how to continuously deform the dimer 
state into the RVB state; as we will show in Sec. lVCl this 
interpolation does correspond to a continuous change of 
the associated local parent Hamiltonian. To this end, we 
define a family of PEPS projectors 



PW = |+)[|0)((02| + (20|) + |1>((12| + (21|)] 
+ ^|-)[|0)((02|-(20|) + |1)((12|-(21|)] 



(12) 



with a two-qubit physical space. Clearly, V{0) = V^ 
and V{1) = 'Pi, both up to local isometrics. By placing 
'P(O) in the tensor network Fig. [21 we therefore obtain 
a smooth interpolation |^(^)) between the RVB and the 
dimer state. Clearly, we can accordingly define P^{0) 
and subsequently 0^{6) = V^{0)V^ ^ which is a contin- 
uous invertible map from rgP^ to rgP^(^). Note that 
while it might seem that the interpolation ([12]) results in 
a discontinuity as 6> ^ 0, the discussion about the equiv- 
alence of RVB and dimer state in the preceding section 
shows that this discontinuity disappears when consider- 
ing whole stars. 

This interpolation can be understood as a way to make 
different dimer configurations more and more orthogonal, 
analogous to what is achieved by decorating the lattice in 
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FIG. 10. a) Lattice used for the numerical calculations. The 
red dots denote P(0) tensors, b) By blocking three spins, we 
transform this to a square lattice, c) Resulting square lattice. 
For OBC, we set the bonds at the boundary to |2) as shown; 
for CBC, we consider a vertical cylinder where left and right 
indices are connected, d) Tensor representation of the scalar 
product (^((9)|^(6>0>. 



Ref. |22|: Overlapping two different dimer configurations 
gives rise to a loop pattern, and the absolute value of their 
overlap is x^~^^, with N the number of loops (including 
length-2 ones) and L their total length (which equals the 
number of lattice sites). Here, x = l/V^ for the RVB and 
X = for the dimer state, and Eq. (p!2]) corresponds to an 



interpolation x{0) = 



1 i-e^ 
V2 1+^2 • 



B. Numerical results 



indices, while the upper and lower outgoing indices are 
still set to |2}. We will denote by Nh and Ny the num- 
ber of sites in each row and column of the square lattice, 
respectively. We will generally consider the case where 
Ny ^ oo for finite N^- In this case, cylindrical bound- 
ary conditions give the same results as periodic boundary 
conditions as long as topological degeneracies are appro- 
priately taken into account, and we will therefore refer 
to them as periodic boundary conditions (PBC) in the 
following. 



1. Correlation functions and gap of the transfer operator 

Firstly, we have studied the decay of correlation func- 
tions along the dimer-RVB interpolation. Correlation 
functions can be computed for lattices of arbitrary size 
using standard methods for the contraction of PEPS J^ 
We have computed various correlation functions, all of 
which show an exponential decay with a slope which 
changes smoothly with 6>, and do not exhibit any sign 
of a phase transition; as an example. Fig. [11] shows the 
connected dimer-dimer correlation function between two 
z <^ z operators [z = diag(l, -1, 1)], where each dimer is 
supported on one square tensor of Fig. [TOb and the two 
squares are in the same row, as a function of the distance 
r between the squares, for different values of 6>, evaluated 
on an infinite system. (This data has been obtained by 
row- wise contraction using infinite Matrix Product States 
and iTEBD, and then evaluating the correlation function 
using the fixed point tensor.) 

The decay of any correlation function is bounded by 
the gap of the transfer operator, i.e., one row in Fig.fTQti. 
We have used exact diagonalization to determine the gap 
of the transfer operator for PBC. The result for N^ - 
4, . . . , 10, as well as the extrapolated data for N^ = oo, is 
shown in Fig. [121 Again, the gap of the transfer operator 
stays finite and converges rapidly for growing Nh- Note 
that the fact that the ground state is a PEPS throughout 



We have numerically studied the interpolation ()T2]) be- 
tween the dimer state (which is a topological Z2 spin liq- 
uid) and the RVB state to determine whether they are in 
the same phase. We use the tensor network representa- 
tion of Fig. m where the interpolating path is character- 
ized by tensors P{0) corresponding to V{0) of Eq. ([T2]) . 
For the numerical calculations, it is particularly conve- 
nient to block three spins, i.e., P{0) tensors, together 
with two E tensors, as indicated by the dashed squares in 
Fig. [TOk : this allows us to rewrite the system on a square 
lattice with three spins per site (Fig. [TOb). We perform 
simulations both for open boundary conditions and for 
cylindrical boundary conditions. For open boundary con- 
ditions (OBC), we set all open indices in the tensor net- 
work to |2} (i.e., there are no outgoing singlets at the 
edges), as indicated in Fig. [TOh . The case of cylindrical 
boundary conditions is obtained by putting the system 
on a vertical cylinder, i.e., connecting the left and right 




FIG. 11. Two-point correlations as a function of the distance 
for different values of 0. The plot shows data for an infinite 
system where in the PEPS contraction,^ the bond dimension 
has been truncated at D = 40 and D = 60, respectively. 
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FIG. 12. Gap of the transfer operator A as a function of 
Q along the RVB-dimer interpolation for Nh = 4, . . . , 10 and 
PBC. As the transfer operator is parity preserving (chang- 
ing) for Nh even (odd), it has two (almost) degenerate 
largest eigenvalues corresponding to different parity sectors 
(for Nh even, there is an exponentially small splitting in 
A/'/j,); the gap A shown is the ratio of the third eigenvalue 
and the average of the two maximal eigenvalues. The data 
for Nh = oo has been obtained by fitting the data for Nh = 
6, . . . , 10 with a + b exp(-cA/'/j,). The inset shows the fit for 
= 0,0.1,0.2,0.3,0.5 (from bottom to top, the curves have 
been shifted vertically in the plot); note that the two right- 
most data points have not been used for the fit. 



the interpolation (and thus satisfies an area law) does not 
preclude critical behavior.— ^2^ 



2. Overlap 



FIG. 13. The overlap change Fn^ , Eq. {[HI), as a function of 
for both open boundary conditions (Nh = 4, 6, 8) and periodic 
boundary conditions (Nh = 4,6,8,10). The solid black and 
red line give the extrapolated data for Nh = oo for open and 
periodic boundaries. The inset shows a zoom of the Nh = oo 
data for < 0.18, where the black (red) crosses correspond to 
open (perodic) boundaries. The details of the scaling analysis 
are given in Fig. 1141 



to arbitrary accuracy for Nh = 4,6,8; this data is then 
used in a finite size scaling analysis to infer f(0) = 
\im Fn^{0). We apply this approach to the OBC 
case. The second approach makes use of the fact that 
lim^^o FNh,N^iO,e) can be expressed as an average over 
correlation functions, as we show in Appendix [Bl Again, 
this average over correlation functions can be computed 
exactly for any value of Ny^ and N^ < 10 (to reach N = 10, 



In order to detect phase transitions which are not re- 
flected in correlation functions, we have also studied the 
rate at which the state changes as we change ^, as quan- 
tified by the overlap (or fidelity). ^^ Concretely, let |^(^)} 
be the state interpolating between |^rvb} = 1^(0)} and 
l^dimer) = 1^(1) }? ^is defined in the preceding section. We 
define the change of overlap 



Fn^^n^O^c) := — 



1 



\{^{0^e)\^{0))f^ 



(^(i9 + e)|^(6> + e))(^(6>)|^((9))_ 

and are interested in the behavior of the limit 

1 



fiO) 



lim lim 



Ar,,Ar^^ooe^o NyNh 



FN.MO.e) (13) 



which describes how quickly the state changes in the ther- 
modynamic limit. 

We have computed this quantity in two different ways. 
In the first approach, we use that for sufficiently small Nh 
we can exactly compute F^^Nvi^^ ^) ^^r arbitrary values 
of A^^ and e by exact row- wise contraction of the tensor 
network Fig. [TOji, and use this to extrapolate the limit 



Fn^{0) = lim lim 



1 



N^^ooe^O NhNy 



FN,,NAe,e) (14) 
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FIG. 14. Finite size scaling analysis for Fn^- The plot shows 
OBC and PBC data for Fn^ ■ Due to the absence of boundary 
effects, the PBC data converges very quickly to the Nh = oo 
value once Nh is beyond the correlation length. We have 
fitted the PBC data for Nh = 6,8, 10 with a + b exp(-cAr^). 
The OBC data, on the other hand, has a 1/Nh leading term 
from the boundary; we have fitted the data for Nh = 4, 6, 8 
with a + b/Nh . Note that the correction for the data points 
for OBC would most likely be concave, suggesting that the 
true Nh = oo value is closer to the PBC data. 
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FIG. 15. Entanglement spectrum of an infinite (vertical) cylinder for decreasing 6, and perimeter Nh = 8. The spectrum 
is shown as a function of the momentum K along the one-dimensional (horizontal) edge. The eigenstates are also labelled 
according to their spin quantum numbers (| representation) according to the symbols in the legend. Panel a) shows the 
integer spin sector of the spectrum, panel b) the half-integer sector; the two sectors are not coupled by the entanglement 
Hamiltonian. Note that for both sectors, the same energy scale has been chosen (i.e., trace-normalizing the boundary states in 
both sectors), corresponding to a local entanglement Hamiltonian.— 



we use that total bond dimension of the transfer opera- 
tor of the P tensor, and thus the overall bond dimension 
in the tensor network in Fig. [TOji, can be reduced from 
D^ = 9 to 6), which again allows us to determine Fn^ 
for Nh = 4,6,8,10, and subsequently apply a finite size 
scaling analysis. We apply this approach to the PBC 
case. 

The resulting data is shown in Fig. [131 together with 
extrapolated data for N^ = oo from the finite size scal- 
ing. (See Fig. [M] for a discussion of the finite size scaling 
analysis.) The extrapolated curves for open and peri- 
odic boundary conditions agree very well. The extrap- 
olated value for f{6) seems to become essentially con- 
stant for < 0.15, which might suggest a non- analytic 
behavior of f(0); however, zooming into this region (in- 
set of Fig. [13]) shows clearly that f{6) decreases again for 
^ < 0.12, showing no evidence for the presence of a phase 
transition. (See Fig. [14] for a discussion of the difference 
observed between the extrapolated OBC and PBC data.) 



3. Entanglement spectrum 

Finally, we have also studied the behavior of the entan- 
glement spectrum along the interpolation for an infinite 
cylinder, using the techniques described in Refs. |2J and 
[25] (in particular, using exact contraction). In Fig. [15] we 
give the entanglement spectra for N^ = 8 for the integer 
and half-integer spin sector (the boundary Hamiltonian 
does not couple the two sectors). While at first, it seems 
that the lowest lying S = ^ and S = 1 states cross with 
the S = ground state at some small 6>, a scaHng anal- 
ysis (Fig. [T6|) suggests that the gap of the entanglement 
Hamiltonian is more likely to vanish at ^ = 0; also, note 



that there is no prior reason to assume that a vanish- 
ing gap in the entanglement spectrum implies any criti- 
cal behavior in the actual system, as the boundary state 
is a thermal state of the entanglement Hamiltonian and 
therefore will have finite correlations even for a critical 
entanglement Hamiltonian.— 

Let us also point out an interesting aspect of the en- 
tanglement spectrum: While the system under consid- 
eration is topological (rigorously provable in an environ- 
ment of the dimer point™), the entanglement Hamilto- 
nian has a unique ground state, i.e., it is gapped, which 
shows that the connection between topologically ordered 
systems and gapless entanglement spectra might be less 
clear than generally believed. Note that the reason that 
the entanglement spectrum is gapped (unlike for the toric 
code, the RK model,— or the mapping of Ref. [20) is due 




FIG. 16. Gap between the lowest S = 1 and S = state 
(left) and S = ^ and S = state (right) of the entanglement 
Hamiltonian for the dimer- RVB interpolation for Nh = 6,8. 
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to the fact that our representation of the dimer state 
involves singlets (i.e., an entangled state of two spins) 
rather than single spins indicating the presence or ab- 
sence of a dimer, and that we split these singlets when 
computing the entanglement spectrum, just as one would 
do for the RVB state itself. 



HAMILTONIANS 



that we have a completely local description of our sys- 
tem; in particular, it is this locality of the description 
which allows us to interpolate between the dimer state 
and the RVB state locally. Beyond that, the particular 
way in which we make the dimer configurations locally 
orthogonal allows us to define the resonance moves Hq 
on a single hexagon (cf. Appendix [C]) . as compared to 
the star required in the usual RK HamiltonianJ^ 



A. Hamiltonian for the dimer state 



Hamiltonian for the RVB state 



As has been proven in Ref.[l6|, every G-injective PEPS 
on a square lattice appears as the ground state of a local 
Hamiltonian, acting on at most 2x2 elementary cells, 
which has a fixed degeneracy determined by the sym- 
metry group. In particular, this result can be directly 
applied to the dimer state to infer that it is the ground 
state of a local Hamiltonian with four-fold degeneracy. 
Alternatively, we can obtain a Hamiltonian for the dimer 
state by using its local isometric equivalence to the toric 
code: The same transformation will also transform the 
toric code Hamiltonian into a local Hamiltonian for the 
dimer state. Note, however, that we need to add extra 
terms to the dimer Hamiltonian which ensure that the 
state is constrained to its local support, as the isometry 
is only defined as a map between the local support of the 
toric code and of the dimer state, respectively. 

In Appendix [Cl we explicitly derive a Hamiltonian for 
the dimer state using the duality with the toric code 
Hamiltonian. It has three types of terms (all of them 
projectors), acting on vertices, triangles, and hexagons, 
respectively: 

• h^ acts on individual vertices. It ensures that there 
is exactly one singlet per vertex. 

• h^ acts inside the triangles (as used in the A-B- 
blocking). It makes sure that each triangle holds 
zero or one singlet. 

• Hq acts on the six vertices adjacent to each 
hexagon. It makes sure that all singlet configu- 
rations around the corresponding star appear with 
equal weight. 

(The exact form of these three terms, together with their 
derivation, can be found in Appendix O) 

How does this Hamiltonian compare to the Rokhsar- 
Kivelson (RK) Hamiltonian for dimer models?—"— The 
resonance terms Hq of our Hamiltonian correspond to 
those of RK Hamiltonian; in both cases, they ensure that 
the ground state is the even weight superposition of all 
dimer configurations. On the other hand, our Hamilto- 
nian has two additional types of terms {h^ and h^) which 
ensure that the system is in a valid dimer configuration; 
these terms are not present in the conventional dimer 
models as they are defined right away on the subspace of 
valid dimer configurations. While this makes our Haml- 
tonian more complicated, it is outweighed by the fact 



The Hamiltonian H^ = Jlhf for the dimer state is 
frustration free^ this is, all Hamiltonian terms — which we 
have chosen to be projectors — annihilate the dimer state, 
^tl^dimer) = 0. We cau make use of this fact to directly 
obtain a corresponding frustration free Hamiltonian for 
|*RVB) = (0^i)®^|*dimer>: Define 



h, ■= oThfoi 



(15) 



where for each i, the product O?^ only involves the stars 
in the covering Fig. [9] which overlap with hf. Then, 
hi is positive semi-definite, and /?^i|^RVB} = 0. More- 
over, for each star k we add a term hk to the Hamilto- 
nian which projects onto the orthogonal complement of 
rg'Px^, and thus restricts the ground state (=zero-energy) 
subspace of the overall Hamiltonian to igV^ . [Note 
that conversely, those terms h\ in H^ which restrict the 
ground state space to (rg'Pi)®^ give /i^ = in Eq. (fT5]) .] 

As O^ is a bijection between rg'Pi and rg'P^, it fol- 
lows that there is a one-to-one correspondence between 
ground states of H^ and oi H ^ Y.hi] m. particular, the 
constructed RVB parent Hamiltonian has a topology- 
dependent degeneracy of the ground space (four-fold on 
the torus), which is obtained from the ground space of 
the dimer Hamiltonian by applying O^ . 

The Hamiltonian constructed this way is rather large: 
In particular, those triangle terms /i^ which lie at the in- 
tersection of three stars give rise to terms acting on three 
stars in Fig.[9l In order to obtain simpler terms, one can 
choose to adapt the framework for parent Hamiltonians 
of G-injective PEPS derived in Ref. 16 to the RVB state, 
which results in a parent Hamiltonian acting on two over- 
lapping stars, i.e., 19 spins; we give the full derivation in 
Appendix [Dl 



C. Hamiltonian for the dimer- RVB interpolation 

The relation between the dimer and RVB Hamilto- 
nian extends to the whole interpolating path introduced 
in Sec. IIV A[ and gives rise to a continuous path of lo- 
cal Hamiltonians interpolating between the dimer model 
and the RVB state: In order to ensure continuity, we 
first replace the dimer Hamiltonian H^ = Xli hf by an 
equivalent Hamiltonian H^ := Xli hf + Y., Uj. Here, hf is 
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the projection of hf onto (rg7^i)'^^*\ where the tensor 
product (S)tv{i) goes over all stars supporting hi and we 
omit vanishing hf^ and 11 j is the projector onto rgT^i 
on star j, where the sum ranges over all stars. As H^ 
is a parent Hamiltonian, the restriction to rgT^i is en- 
sured by the terms in H^ locally; thus, H^ has the same 
ground state space as H^ with a spectal gap above, 
and we can smoothly interpolate between H^ and H^ 
without changing the ground state space or closing the 
gap. For H^^ it is now straightforward to construct a 
continuous interpolating path H{0) = Y.ihi{0) + llj{0)^ 
where hi{e) = O^ {e)hiO^ {0) , and n^(l9) projects onto 
rg P^(6>) for star j. 



VI. CONCLUSIONS 

In this paper, we have applied the PEPS formalism 
to the study of the Resonating Valence Bond states and 
dimer models. We have discussed PEPS representations 
of the RVB and the dimer model and studied their struc- 
ture and relation. In particular, we have given a local 
unitary mapping between the dimer model and the toric 
code; furthermore, by defining the dimer state with lo- 
cally orthogonal dimers, we were able to devise a local re- 
versible mapping between the dimer state and the RVB 
state which allowed us to prove that the RVB state is 
the four-fold degenerate ground state of a local parent 
Hamiltonian for any finite lattice. Subsequently, this al- 
lowed us to devise a smooth interpolation between the 
dimer state and the RVB state and the corresponding 
Hamiltonians. We have studied this interpolating path 
numerically, considering the behavior of correlation func- 
tions, the rate at which the ground state changes, and the 
entanglement spectrum, and have found that all of these 
quantities behave smoothly and show no sign of a phase 
transition. 

Our results make heavy use of the formalism of PEPS 
and their associated parent Hamiltonians, and in particu- 
lar of G-injectivity and G-isometry, which we generalize 
in order to assess the RVB state on the kagome lattice. 



Similarly, the PEPS representation of the dimer-RVB 
interpolation allowed for efficient numerical simulations, 
enabling us to study the phase of the RVB state. We 
believe that these techniques will be of further use in the 
study of related systems. In particular, the results can 
be generalized to other lattices: Firstly, the PEPS de- 
scription of RVB states applies to arbitrary graphs. All 
lattices with the "linear independence property" ^i^^i^^ 
are G-injective, and whenever the linearly independent 
blocks allow cover the lattice up to disconnected patches 
(such as in Fig. [9]), this allows to interpolate between the 
RVB and the corresponding dimer state. (This is the 
case, for instance, for the square-octagon or the star lat- 
tice in Ref.— , whereas for the hexagonal lattice it appears 
that no such covering can be found.) Further, if the dimer 
model can be expressed in terms of Z2-injective tensors, 
this allows to conclude that it is equivalent to the toric 
code. (This is the case e.g. for the star lattice, but not for 
the square-octagon lattice.) Also note that our findings 
are not restricted to spin-| SU(2) singlets as dimers, but 
can be generalized to higher dimensional singlets or any 
other state, as long as Z2-injectivity can be established. 
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Appendix A: Different PEPS representations for the 
RVB and dimer state 

In this appendix, we give an overview of different PEPS 
representations of the RVB and dimer state and how they 
are related. We wih start from the representation intro- 
duced in Ref 17, and subsequently show how to derive 
from it the representation used in this paper. 



I 0]. The two-dimensional subspace, (5 = span{|0}, |1}}, 
provides the spin-^ degree of freedom which holds the 
singlets |01}-|10} used to construct the RVB state, while 
the third level, |22}, is used as a tag indicating that there 
is no singlet along that edge, synchronizing this choice for 
both ends. In order to understand what happens when 
we apply Va to the bonds |cj}, let us rewrite Va as 



Va = Va V. 



4,1 



1. PEPS representation of the RVB 

We first explain the PEPS representation of the RVB 
state introduced in Ref.llTI, which is illustrated in Fig.fTTl 
We place states 



a; 



|01)-|10) + |22)gC^(8)( 



(Al) 



along all edges of the lattice (observing its orientation); 
this associates four three-level systems (qutrits) with 
each vertex. Then, we apply the following map to the 
qutrits at each vertex, which maps them to one physical 
qubit: 



P4 



=|0)((0222| + (2022| + (2202| + (2220|) 
+ |1)((1222| + (2122| + (2212| + (2221|) 



(A2) 



= Z(|o)(ok 

fe=l 



|l)(l|fc)0(222|/fc. 



Here, in the second formulation the sum runs over the 
four virtual systems /c, (0|/c and (l|/c act on the virtual 
system /c, and (222|//g acts on all virtual systems but k. 

What is the intuition underlying this construction? 
The Hilbert space holding the bond state, \uj) = |01} - 
1 10} + |22}, Eq. (|A1|) . should be understood as a direct 
sum C^©C^ [with a corresponding SU(2) representation 




FIG. 17. PEPS representation for the RVB. Place maximally 
entangled "bonds" |01) - |10) + |22) along the edges of the lat- 
tice (blue), and subsequently apply the linear map ("PEPS 
projector") V4 (red circle), Eq. ()A2|) . which maps the four 
qutrits at each vertex (encircled) to a spin-| degree of free- 
dom. Alternatively, we can obtain the (orthogonal) dimer 
state by replacing V4 with 7^4,1, Eq. ()A3|) : 7^4,1 ensures that 
only one qutrit per vertex holds a singlet, but keeps the full 
Hilbert space, thus ensuring local orthogonality of different 
dimer configurations. 



where V4,^i is a projector acting within the virtual space 
as 

V4,,= t^e,k^\222)/k{222\/k, (A3) 

k=l 

with te,k the projector onto the subspace 6 on site k. 
Thus, 7^4,1 projects one virtual qutrit onto the "singlet 
space" (3, while the others are projected onto the |2) sub- 
space. Due to the form of |a;), it follows that both sides of 
each bond have to be projected onto the same subspace, 
i.e., either the singlet is selected, or both ends of the bond 
need to be in the |2) state. From these two observations, 
it follows that by applying V4^_l to the bonds |cj), we ob- 
tain configurations with singlets along the edges, such 
that exactly one singlet is incident to each vertex, while 
all other edges are in the "unused" state |22). This way, 
all possible dimer configurations are orthogonal — we can 
infer the location of the dimers by measuring whether we 
are in the 6 subspace or in the |2) state. Moreover, the 
final state has an equal weight of all dimer configurations, 
as |01)-|10) and |22) appear with the same relative weight 
everywhere. Thus, it follows that the state obtained by 
applying the maps V^^i to the bonds \uj) is a realization 
of the (orthogonal) dimer state. 

It is now straightforward to see that applying V4 to 
this dimer state gives the actual (non-orthogonal) RVB 
state: V4 keeps the singlet degree of freedom at every 
vertex, but it (coherently) erases the information about 
the edge the singlet is associated to; as the phases be- 
tween different dimer locations are chosen to be +1, we 
indeed obtain an equal weight superposition of all dimer 
coverings with singlets, i.e., the RVB state. 



2. Simplified PEPS representation 

We will now prove that the PEPS representation of the 
RVB state Eqs. ()AllA2p is equivalent to the form intro- 
duced in Sec. IIIBi We start by rewriting V4,, Eq. (jAip . 
as 



V4 = V{Q^Q) . 



(A4) 



Here, V • C^ <S) C^ ^ C'^ is defined just as 7^4, but acting 
only on two virtual systems, 

P = |0)((02| + (20|) + |l)((12| + (21|) 
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FIG. 18. Simplified PEPS construction for the RVB and 
dimer state, a) Start from the representation Fig. [T7] and 
rewrite P4 in two layers, Va = V(Q <S> Q), Eq. (lA4l) . The 
Q act as indicated by the green circles, and make sure at 
most one singlet degree of freedom is kept; they thus return 
one qutrit each. V2 makes sure exactly one singlet is kept, 
which can come from either of the two Q's. b) The simpli- 
fied PEPS construction: Applying three Q's to three bonds 
around a loop results in a state |e) which holds at most one 
singlet, Eq. (|A5|) . and which is placed across every triangle. 
Subsequently, V is applied to pick one of the singlets at each 
vertex. 



C^ -^ C^ is defined as 
|1)((12| + (21|) + |2)(22|; 



[cf. Eq. (j2j)], and Q : C^ ^ 
QH0)((02| + (20|) 



it is straightforward to check that (|A4p holds. We choose 
to apply the Q's as depicted in Fig. [T8k : Then, applying 
Q®^ to \uj)®^ across each triangle yields a 3-qutrit state 



k)= E e^jk\ijk) + \222) 

i,j,k=0 



(A5) 



where eijk is the completely antisymmetric tensor with 
5oi2 = 1, where i, j, and k are oriented clockwise. Thus, 
we obtain a PEPS construction for the RVB state where 
we place tripartite bond states \e) across triangles, and 
apply V to each adjacent pair of qutrits, as shown in 
Fig. UHlo. 

As already explained in Sec. Ill CI this PEPS repre- 
sentation of the RVB state can again be understood as 
arising from a PEPS representation of the dimer state, 
V = V V^ with V^ ^ t& ® |2)(2| + |2)(2| le- Note 
that different dimer coverings are again locally orthog- 
onal, although in a different fashion — they become or- 
thogonal on triangles. Note that this dimer state can 
be obtained from the one discussed in the previous sec- 
tion [AT] using that V^{Q® Q) = (7^ 7^)7^4,l, where 
7^ = 1© (8) (2| + (2| (8) le + |2)(22|: While 7^ removes some 
degrees of freedom, different dimer configurations remain 
orthogonal on triangles. 



Appendix B: Change of overlap and correlation 
functions 



In this appendix, we show how to compute the change 
of the overlap of the state along the dimer-RVB inter- 
polation from two-point correlation functions; it can be 



seen as a Hamiltonian-free PEPS version of the result 
that the change of the overlap can be expressed as an in- 
tegral over imaginary time correlation functions. ^^ While 
we will illustrate the derivation for the particular case 
of the dimer-RVB interpolation, it will equally apply to 
most other PEPS interpolations. 
The interpolating path is given by 



\^9)=pr\^), 



with V0 as in Eq. ((T2|) . and where \D,) = l^fdimer) 
dimer state. Ve can be spht as 

Pe = P^+ OP- , 



(Bl) 
is the 



(B2) 



where P+ and P_ are projectors independent of ^, P+ + 
P_ = 1; the following derivation holds for any interpola- 
tion of the form Eqs. (|BllB2p . Note that from (|B2]), we 
have that 

Ve+de = {t + lP_dO)V9 . (B3) 

We want to compute the change of the normalized overlap 



0(0,0 + dO) 



\/{^e\^e)\/{^e+de\^e+de) 



As the first-order change to O{0,d0) can be at most a 
(unphysical) phase change, we are interested in the sec- 
ond order in dO. We write \ipe+de) - \^e)^\di^e)-, and from 
Eqs. (|Bip and (|B3p , it follows that \dilJo) can be expanded 
in orders of dO. 

For notational convenience, let us define 

B = {d^0\di>e)/M . 

Note that since both \Q) and P0 in (jBip are real, iV'e) is 
real, and thus A is real; we will use this in the following 
derivation. Using the above definitions (where A is first 
order and B second order in dO) we have (neglecting third 
and higher order terms) 



0(6,6 + dO) 



^fil + A) 



^f^l + i2A + B) 
(l + A) 



2A + B 3(2A)2 
1 + — ^^ — — 



= i-Ub-a^) . 



Clearly, 



dPe, 



l#^> = E ^p!'''-'>\n)d0 = Z -^P-Mde , 

where the sum is over all possible positions i of the hPl- 
(As A only appears as A^ in the overlap, we can neglect 
the second order term.) It follows that 



^^i^{MPi\i>0 



1 



A/" 



d6^-Z{P-hde 



14 



where we have used the notation 
On the other hand, to second order 

ij 

and thus, the second-order term in the expansion 

0(e,e + d0) = 1 - IgeNidef 

can be expressed as an average over connected correlation 
functions: 



a) 



B-A^ 1 



NidO)"^ O'^N 



Y.{P-P'- 



(B4) 



where P! '= Pl-{Pl)e^ and N is the number of sites. This 
derivation shows that for PEPS interpolations, we can 
infer the change of the overlap from two-point correlators; 
this argument can be extended to relate the p-th order 
expansion coefficient of the overlap to p-point correlators. 

For systems with finite correlation length ^, the sum 
in Eq. (|B4|) is proportional to ^'^N, and thus, ge - ^^/6>^. 
Therefore, the only way ge can diverge when the corre- 
lation functions decay exponentially is if ^ ^ 0. Note, 
however, that for 6 = 0, {PlPl)e = 0, so whether ge di- 
verges at ^ ^ depends on the rate at which the latter 
vanishes. 

Note that Eq. (|B4[) does not apply when ^ = 0. For 
that case, one can easily check that {ilje=o\dilJe=o) = and 
thus A = 0. On the other hand, {dilJe=o\dtlJe=o) / {dOy is 
the norm of the PEPS where Ve=o = P+ has been replaced 
by P_ at a single site (summed over all sites), which im- 
mediately shows how to compute go = B/N{dO)^. 



Appendix C: A simple parent Hamiltonian for the 
orthogonal RVB 

In this appendix, we show how to directly derive a 
parent Hamiltonian for the dimer state (in the represen- 
tation of Sec. IIICp by relating it to the toric code on the 
same lattice. 

Let us start by considering the toric code model on a 
kagome lattice, where the spins sit on the edges of the 
lattice, as illustrated in Fig. [T9k . The toric code Hamil- 
tonian has two types of terms: First, it has terms acting 
on the four qubits across each vertex, as 

and second, it has terms acting across plaquettes — 
triangles and hexagons — as 



/-I 
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FIG. 19. a) Tensor network and Hamiltonian for the toric 
code on the kagome lattice. Two types of terms ensure the 
parity constraints over the triangular and hexagonal plaque- 
ttes, while the vertex term enforces an equal weight superpo- 
sition of all even parity configurations. Note that the -Z"®^ 
acting on the triangles can be understood as enforcing each 
individual tensor to be in the correct subspace. b) Hamilto- 
nian terms obtained by transforming each tensor unitarily into 
a projector onto the symmetric subspace. Z®^ still enforces 
that the tensor is in the correct subspace; further, there is a 
Z®^ at each vertex enforcing equality of the adjacent qubits — 
translated to the dimer state, this will enforce that there is 
exactly one singlet adjacent to each vertex. Finally, the X®^^ 
term around the hexagon makes sure that all such config- 
urations appear with equal weight — it will make the dimer 
configurations resonate in the orthogonal RVB state. (Note 
that this model can also be understood as a toric code model 
on a "decorated" hexagonal lattice with two sites per edge 
and a vertex between, and X stabilizers around plaquettes 
and Z stabilizers across vertices.) 



where X and Z are the Pauli matrices. For future pur- 
poses, we have introduced different labels for the two 
types of triangles: h^^ refers to an upward pointing ten- 
sor (sketched blue in Fig.fTOk). and thus to a downwards 
triangle, and vice versa. In analogy to the square lattice 
toric code, Eq. (|6|), we can write the toric code state with 
triangular Z2-invariant tensors 



V' 



q;,/3,7 



The states {fj\ on the virtual level are again in the x 
eigenbasis, which makes them invariant under Z®^. The 
tensors are marked blue in Fig. [T9k . with the virtual in- 
dices associated to vertices of the kagome lattice. V^^ 
can be transformed by a unitary acting on the physical 
system to the projector onto the invariant subspace, 

where the physical qubits are now associated to vertices, 
cf. Fig. [T9b . It can be easily checked that the action 
of the Hamiltonian terms in the transformed basis is as 
follows (illustrated in Fig. [T9b): 

enforces that the two qubits per vertex are in the same 
state (namely equal to the value of the virtual index); 



/^r 



hi"" = -z^ 



TC 



= -z^ 
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acts on the 12 qubits around each hexagon and ensures 
that all bond configuration appear with the same prob- 
ability; finally, in the original representation, h^^ = h^ 
enforced that the qubits at each tensor are in the proper 
subspace, it thus has to be replaced by the corresponding 
term for the transformed tensor. 

The three terms are depicted in Fig. [T9b . 

Let us now consider what happens if we change all A- 
type tensors (upwards pointing triangles) to projectors 
onto the antisymmetric subspace, 

while the 5-type tensors (downwards pointing triangles) 
remain of the form 

Clearly, this gives rise to a different triangular Hamilto- 
nian for all A-type tensors, 

Tdef _ ry®?, 




which enforces odd parity, while we keep N^^ = h^^ for 
the B-type tensors. All other Hamiltonian terms remain 
unchanged, since we still need to ensure that the qubits 
on both sides of the vertex are equal {h'^^^ = h^^), and 
that all configurations appear in an equal weight super- 
position {h"^^ =Kly^). 

The projectors V"^^ and V^^^ can now be transformed 
to the A and B tensors for the dimer state by local isome- 
trics, which allows us to explicitly construct its parent 
Hamiltonian: First, the terms h^^ and h^^ enfore that 
the sites at each tensor are in the correct subspace; they 
are thus mapped to a local term h^ = Ha which have the 
four states of Fig. [71 as ground states. The vertex term 
h"^^^ is transformed to 

/i, = -(le0|2)(2| + |2)(2|0le) , 

(5 = span {|0}, |1}}, which ensures there is only one singlet 
adjacent to each vertex. Finally, Jiq^ is transformed to 
a term which swaps the singlet/no-singlet configuration 
for the 12 qutrits adjacent to each hexagon, observing 
the orientation of the singlets: This is achieved by an 
operator 

where E acts on the two qutrits in each triangle adjacent 
to the hexagon as 

E = |22)((01| - (10|) - |20)(02| - |21)(12| + h.c. , (CI) 

with proper orientation of the singlet — the first term 
changes between a singlet and no singlet along the edge, 
and the latter two terms fiip the location of a (outwards- 
pointing) singlet in the triangle; the negative sign is due 
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FIG. 20. Resonance moves in the dimer state, a) Resonance 
move implemented by a single S, Eq. (|Cip : S acts on the 
lower two qutrits only and flips the two configurations sin- 
glet/no singlet, and the two configurations singlet left /singlet 
right, a) Applying S®^ around the hexagon gives a resonance 
between pairs of dimer configurations around the hexagon 
(marked red and green, respectively). Note that any such 
pair of dimer configuration can be described by a closed loop 
of even length on the star, cf. Ref. |8|. 



to the fact that the two singlets are oriented differently. 
The action of S is illustrated in Fig. [2Qk . Note that 
this S®^, although it acts on a hexagon only, in fact 
achieves resonance between all loop configurations on the 
star built around the hexagon, as illustrated in Fig. [2Qb 
for a particular instance. 



Appendix D: C— inject ive toolbox and parent 
Hamiltonian for the RVB state 



In the main text, we have seen how to manipulate the 
RVB state in a kagome lattice, to connect it with the 
dimer state and then with the toric code in such a way 
that we can apply the results about Z2-injective PEPS 
given in Ref. [l6| to give a local frustration free Hamil- 
tonian with topological features for the RVB state. An- 
other possible approach is to modify the toolbox for G- 
injective PEPS developed in Ref. [l6| to adapt it to lat- 
tices different than the square one, even frustrated ones. 
Since this route is necessarily technical, and assumes the 
reader is familiar with Ref. |l6|, we only sketch it in this 
appendix. The payoff: a more local Hamiltonian — acting 
on two partially overlapping stars, that is, 19 spins — and 
the possibility of applying the G-injectivity toolbox to 
rather general lattices in the future. 

For simplicity, we will stick to Z2-injectivity and the 
kagome lattice, but it is straighforward how to extend the 
ideas to general finite groups and general lattices. We 
will work with the tensor network picture of the RVB 
state given in Fig. IDd, where the RVB state is recov- 
ered from contracting two types of tensors P and E. 
In the first part of this appendix. Sections ID II to ID3( 
we will generalize the three key ingredients of Ref. [l6|: 
(i) showing that Z2-injectivity is stable while growing 
the lattice, (ii) showing the "intersection property" for 
good-shaped Hamiltonians — as the one acting on three 
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FIG. 21. Graphical proof that Z2-injectivity is stable as grow- 
ing a region. We illustrate it in the case of growing from one 
to two partially overlapping stars. Our aim is to prove the ex- 
istence of a left-inverse on the invariant subspace on the two 
stars, using the existence of the left-inverse on a single star. 
Since we are considering the tensor network representation of 
the RVB state made out of tensors P — with physical index — 
and triangular links E without physical index, given a region, 
all one needs is to clarify whether the outgoing auxiliary in- 
dices come from a P tensor or from an E tensor. In our case, 
we start with two overlapping stars with outgoing indices of 
type P and divide them in three regions. The middle one 
(bow-tie) with all outgoing indices of type P and the rest. 
By gathering indices together, we have figure (a), where the 
circle represents the bow-tie and each square one of the other 
two regions. The crucial property is that the circle together 
with any of the two squares is Z2-injective. We proceed as 
follows. We apply the inverse to the circle and left square, 
getting Figure (b), where Z now means ®Z. The next step, 
illustrated in (c) is to build a new circle tensor attached to 
the right square tensor, and passing the Z^ through using the 
symmetry of the circle tensor. Now we invert the remaining 
two tensors getting (d). By tracing the middle (extra) line 
we obtain 5rs and therefore finish the argument since we have 
constructed the desired inverse. 





zizlzl 
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FIG. 22. Graphical proof of the intersection property for the 
case hi + h^ where hi acts on the 3 stars colored in blue-red 
in (a) and /i2 acts on those colored in red-green. It is easy 
to see that, by separating the common part (2 red stars) and 
the rest, with the convention that the common part has all 
outgoing indices of type P, ker hi ® Tigreen is given by figure 
(b), where N is an arbitrary tensor, whereas T^biue (8)ker/i2 
is given by figure (c), where M is arbitrary. Here Tigreen and 
'Hbiue represent the Hilbert spaces associated to the physical 
indices in the green and blue regions, respectively. As before, 
we have gathered together all virtual (as well as the physical) 
indices within each type of tensor. Given a state |'0) in the 
intersection of both kernels, there exist N and M such that 
\ip) can be expressed both in the form (b) and (c). In the 
equality (b) = (c) we invert the circle tensor, obtaining the 
equality given in figure (d). After that, we apply an auxiliary 
extra circle tensor, then the inverse of the red-green tensors (3 
stars), which is also Z2-injective by the previous section, and 
then trace the extra virtual system. Since one can assume 
wlog that both M and N have the symmetry (g)Z, we can also 
eliminate the sums. After all that we arrive at figure (e) for 
some tensor L. Substituting (e) in (c) gives the desired form 
for \ip) and hence the result. 



overlapping stars — and (iii) showing that we recover the 
desired ground space while closing the boundaries. All 
that will produce a local frustration free Hamiltonian for 
the RVB state with local terms acting on three overlap- 
ping stars and 4- fold degeneracy. Finally, in Section ID4| 
we will see how to reduce the size of the local Hamilto- 
nian to two overlapping stars. For that, we will need the 
concept of linear independence introduced by Seidel in 
Ref. i. 



1. Growing Z2-injectivity 

We recall from the main text that a region in a tensor 
network is Z2-injective if it is Z2-invariant, as illustrated 
in Fig. [5l and in addition one can left-invert the tensor 
in that region (understood as a map from the auxiliary 



to the physical system) to get the projection onto the in- 
variant subspace. In Figure [8] of the main text, the basic 
Z2-injective block — one star — is depicted. In Figure [21] 
we illustrate how to grow the Z2-injective property from 
one to two partially overlapping stars, but the procedure 
can be applied to grow arbitrarily large regions. The 
only basic property we need for that is the possibility of 
decomposing the final tensor as in Figure [2Tb . and that 
this property is preserved as we grow the region. The 
general result is that for regions that are grown keeping 
this property (it is trivial that one can grow in this way 
all regions that are required to be Z2-injective later on 
in this appendix), there is an inverse giving the projector 
onto the invariant subspace: 



1 1 
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where |t| is the number of tensors, both of P and E type, 
|6| is the number of outgoing indices in the region, and 



Z:=\ 



-1 



As commented in the main text, the extra sign one ob- 
tains depending on the parity of the number of elemen- 
tary tensors in the region only affects the definition of 
"invariant subspace", which can be the symmetric {\t\ 
even) or the antisymmetric one {\t\ odd). 



2. Intersection property 

We recall from the main text and Ref. [16! that given 
a PEPS, one can construct parent Hamiltonians H = 
Ei? hn where the local term hn only acts on a region 
R. The way to construct them is simple. Each Hr is the 
orthogonal projector (or any positive semidefinite opera- 
tor) whose kernel coincides with the range of the PEPS 
tensor on that region (seen as a map from the auxiliary 
to the physical indices in the region). This Hamiltonian 
is > 0, frustration free and has the initial PEPS as ground 
state. In Ref. [I^ we noticed how in the square lattice, 
G-injectivity preserves the structure of the parent Hamil- 
tonian as we grow it properly (i.e., by translational in- 
variance) from a single local operator. That is, the kernel 
of Ei?c7^ hn is exactly the range of the PEPS tensor in 
the larger region IZ. This property is known as the "in- 
tersection property" since it is related to the way the 
kernels of Hr for different regions intersect. We will il- 
lustrate here the intersection property for the case /ii + /12 




FIG. 23. The kagome lattice with periodic boundary condi- 
tions can be decomposed as in (b), where all outgoing indices 
on the black regions are of type P, whereas all outgoing in- 
dices of the blue regions are of type E. All tensor in (b) have 
an additional physical index "going out of the paper". A pos- 
sible way to obtain such a decomposition is illustrated in (a). 
The black regions are Z2-injective, the same happens if one 
takes two black regions with the intermediate red or green 
one or if one takes the four black regions with the interme- 
diate red, green and blue ones. The blue region is the one 
corresponding to the cross points. Note that the black and 
the blue only interact through a red or green region, exactly 
as required in (b). 



where hi acts on the 3 stars colored in blue-red in Figure 
[22k and /12 acts on those colored in red-green. The only 
property one needs to make the argument work is that 
the common interaction area is Z2-injective and that the 
non-common areas are not directly connected. The proof 
is graphically given in Fig. [22l 



3. Closing the boundaries 

It is not difficult to see that, if one considers the RVB 
state as a PEPS in a torus of arbitrary size, by block- 
ing tensors and gathering indices together, one has the 
tensor network given in Figure [23b. It may seem at first 
sight that we are back to a square-lattice PEPS. How- 
ever, there is a crucial difference. Not all tensors in the 
square decomposition of Figure [23b are Z2-injective. We 
do require that the four overlapping regions T^i , . . . , IZ^ 
resulting from removing one row and one column with 
an intersecting blue tensor contain all possible 3-star re- 
gions and that they can be constructed by growing 3-star 
regions keeping the intersection property commented on 
last section. Figure [23k illustrates a possible way of get- 
ting that in the surrounding of a blue tensor. We will 
now use the trivial fact that, given the frustration free 
Hamiltonian 



H : 



E 



(Dl) 



se{3-star} 



resulting from summing the local projectors correspond- 
ing to each 3-star region in the torus, its ground space 



coincides with n'^^-^kei H^^ 



where Hk 



^sc7^/ 
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FIG. 24. Any state which is in the ground state of the parent 
Hamiltonian when we close the boundaries in a periodic way 
(where we close the boundaries as indicated in Figure [23h ) 
can be represented graphically in any of the forms given in 
the figure. As in Figure [23l each of the tensors, including the 
"boundary conditions" have a physical index going "out of 
the paper" . 
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Wr,s) = 




FIG. 25. The ground space of the Hamiltonian with periodic 
boundary conditions has dimension four and is spanned by 
the states |'0r,s) in the figure. As above, Z means ®Z. 



Hamiltonian corresponding to the region IZk- By the in- 
tersection property proved in last section, the associated 
kernels are given by Figure [24l 

The argument now is a modification of Theorem 5.5 of 
Ref.[l6|. We include a sketch for the sake of completeness. 
The aim is to show that the ground space of (|D1[) is 
spanned by the four states \il^r,s)i r^s = 0,1, depicted in 
Figure [251 

Using the two horizontal equalities of Figure [24] and 
reasoning as in the previous sections (inverting the two 
columns which are Z2-injective, growing new tensors 
around the remaining column, inverting them and taking 
traces), we obtain that the state is of the form of Figure 
[26k (from the upper equality in Fig. [24]) . as well as of the 
form of Figure [26b (from the lower equality). The verti- 
cal equalities then give that Fig. [26k equals Fig. [26b . 

We can now invert the two rows containing black ten- 
sors. The way to do it is to invert the two black tensors 
with one intermediate red one, then construct two new 





FIG. 27. The figure illustrates the first step in the reduction 
to show that the parent Hamiltonian on two overlapping stars 
suffices to give the RVB state as a unique ground state - up to 
the string operators given in Figure [25] which induce a 4- fold 
degeneracy. The left (right) upper part of (b) represents the 
kernel of the Hamiltonian grown from 2-star regions which 
acts on the stars with a dot (cross) . The circles represent the 
stars with the same color in (a), with all outgoing indices of 
type P. The squares represent the rest of (a), having outgoing 
indices both of type P (those connecting with the boundary) 
and type E (those connecting with the circles) . As in the pre- 
vious figures, all bonds connecting the same objects have been 
gathered together in the same bond. Also, all tensors have 
an additional physical index which "goes out of the paper" 
The boundary conditions in the upper part of (b) have also a 
physical index corresponding to the part of (a) with the same 
color. The intersection of these two kernels, which is the ker- 
nel of He = Y,x o ^2, is the sum of the two subspaces depicted 
in the lower part of (b). Needless to say that the boundaries 
of the lower part of (b) do not have physical index. 



black tensors around the remaining red one, invert again 
and take traces. This even allows us to select a particu- 
lar r = ro, this way removing the sum over r. We end up 
in Figure [26h and from there, with similar techniques we 
show that the given state must be a linear combination 
of states of the form given in Figure [25] Since clearly 
ah such states belong to the ground space of the Hamil- 
tonian, and all of them are linearly independent, they 
form a basis of the ground space. Note that, as in the 
toric code, all four ground states are related by non-local 
string operators on the virtual level. Indeed, the action 
of these virtual operators is nothing but changing the 
sign (or flipping the direction) of the dimers crossing the 
links in which they lie. Note also that the same proce- 
dure and conclusions can be obtained starting from the 
parent Hamiltonian of a basic region different from the 
3-star one, as long as one can grow properly the lattice 
starting from this region while keeping the intersection 
property. Another example of such a basic region is the 
7-star region of Figure [27k . 



A 2-star Hamiltonian 



FIG. 26. Intermediate step in the proof that the ground space 
of the Hamiltonian is spanned by the vectors given in Figure 
[25] above. 



The first step to show that a 2-star Hamiltonian suffices 
is sketched in Figure [27] Starting from the parent Hamil- 
tonian for 2-star regions, and using the intersection prop- 
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FIG. 28. The figure illustrates the second step in the re- 
duction to a 2-star Hamiltonian. (a) represents the space S, 
coming from the lower-right part of Figure [27b . each state 
of which we want to penalize by the 2-star term marked by 
squares in Figure [27k. All outgoing indices are of type P. (b) 
illustrates the idea of "completing the path": for each set of 
free sites -marked by a circle-, there exist 2 disjoint dimer 
coverings (in green and violet respectively) compatible with 
the choice of free sites, which in addition form a Hamiltonian- 
cycle in the remaining graph. 



erty, one gets that the Hamiltonian acting on the four 
stars marked with a dot (cross) has a kernel equal to the 
upper-left (upper-right) subspace in Fig. [27b . Therefore, 
the kernel of the Hamiltonian Hq = E^ o ^2 — the sum 
of all terms acting on two overlapping stars marked by 
both circles or crosses — is the intersection of those sub- 
spaces. Take a vector in the intersection. Then we can 
reason with the very same techniques as above. First we 
invert the circles, then grow again circles at both sides of 
the blue square and invert the joint tensor circle-square- 
circle. Closing the unnecessary indices by taking traces 
we arrive at the fact that the original tensor indeed be- 
longs to the sum of the subspaces in the lower part of 
Figure [27b . Since clearly each one of these subspaces is 
contained in the kernel of Hq^ the lower part of Figure 
[27l is indeed equal to the kernel of Hq . 

The second step is to add a 2-star term including the 
middle star, for instance, the one marked by squares in 
Figure [27b. We will show that this term penalizes all the 
states in the lower-right subspace of Figure [27b . That 
is, naj^o,nkeri7^2 = ker{X!a^ o n ^2} equals the lower-left 
subspace in Figure [27b. where the intersection/sum runs 
on all 2-star Hamiltonians in the 7-star region given in 
Figure [27b . This procedure illustrates how from 2-star 



Hamiltonians one reaches the desired ground space in the 
7-star region of Figure [27b. Since this ground space is the 
same one obtains growing from 3-stars, we are done. 

Let us first explain the intuition behind the following 
reasoning. For this, we will use the PEPS representa- 
tion of the RVB of Fig. [HI The intuition is that the Z 
appearing in the lower-right part of Figure [27b — which 
denotes nothing but a Z®^ in the 4 bonds connecting the 
corresponding regions — cannot escape the middle star: 
That is, no matter how we use the gauge symmetry of 
the tensors in the PEPS, there will be one, and only one, 
Z in one of the bonds of the middle hexagon. This "exci- 
tation" can be then detected by a 2-star term including 
the middle star. To make it rigorous we will need the 
linear-independence property, proven by Seidel in Ref. S 
The sum ^oiWiD)) <S)l-LD,hee} is direct. Here, D means 
a fixed dimer covering of the 2-stars, \a{D)) the asso- 
ciated tensor product of singlets, and HD,hee the total 
Hilbert space of those free sites which do not have any 
dimer. The sum runs on all possible dimer coverings of 
the 2-stars. 

Let us now take the Hamiltonian acting on the stars 
marked with a square in Figure [27b . Using the Z2- 
injectivity, its kernel is exactly given hy JC = ^bi^-b ^ 
Evaiid k(^)}}5 where 6 is a fixed choice of free sites — 
with an odd number of terms to allow dimer coverings — 
and the sum runs over all possible dimer coverings having 
b as the set of free sites. We are then finished if we can 
show that the subspace S given by Figure [28b has trivial 
intersection with /C, or equivalent ly, that the sum tS /C 
is direct. By the very definition of the PEPS tensor, S 
is given exactly by 5 = ®b{Hb ® Evaiid ^k(^))} where 
Z changes the sign of \cf{D)) when D has one dimer on 
a link marked by Z in Figure [28b and only in this case. 
By the linear independence property of Seidel, the sum 
5 H will be direct as long as we can guarantee that for 
each 6, in Y.va\\d^W{^)) there will be terms in the sum 
with dimers on the Z-links and terms without dimers on 
the Z-links. This last statement can be easily checked by 
"completing the path", as illustrated in Figure [28b. The 
fact that for each 6, there exists a Hamiltonian-cycle of 
alternating colors for the sites not in h can be easily seen 
as follows. We can have paths of alternating colors in 
the desired vertices in all but the middle bow-tie. Since 
h contains an odd number of sites, at one side of the bow 
tie we arrive with the same color whereas at the other 
side we arrive with different colors. Exactly as in Figure 
we can complete the cycle within the bow-tie. 
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